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Abstract. We study a natural online variant of the replacement path 
problem. The replacement path problem asks to find for a given graph 
G — (V,E), two designated vertices s,t £ V and a shortest s-t path 
P in G, a replacement path P e for every edge e on the path P. The 
replacement path P e is simply a shortest s-t path in the graph, which 
avoids the failed edge e. We adapt this problem to deal with the natural 
scenario, that the edge which failed is not known at the time of solution 
implementation. Instead, our problem assumes that the identity of the 
failed edge only becomes available when the routing mechanism tries to 
cross the edge. This situation is motivated by applications in distributed 
networks, where information about recent changes in the network is only 
stored locally, and fault-tolerant optimization, where an adversary tries 
to delay the discovery of the materialized scenario as much as possible. 
Consequently, we define the online replacement path problem, which asks 
to find a nominal s-t path Q and detours Q e for every edge on the path Q, 
such that the worst-case arrival time at the destination is minimized. Our 
main contribution is a label setting algorithm, which solves the problem 
in undirected graphs in time 0(m log n) and linear space for all sources 
and a single destination. We also present algorithms for extensions of the 
model to any bounded number of failed edges. 



1 Introduction 



Modeling the effects of limited reliability of networks in modern routing schemes 
is important from the point of view of most applications. It is often unrealistic to 
assume that the nominal network known at the stage of decision making will be 
available in its entirety at the stage of solution implementation. Several research 
directions have emerged as a result. The main paradigm in most works is to ob- 
tain a certain 'fault-tolerant' or 'redundant' solution, which takes into account 
a certain set of likely network realizations at the implementation phase. One 
important example is the replacement path problem (RP) [18111] . The input in 
RP is a nominal network given as a graph G — (V, E), a source s, a destination 



t and one shortest s-t path P. The goal is to find for every edge e on the path 
P, a shortest path in G, which does not use the edge e. RP attempts to model 
the situation in which any link in the network may fail before the routing pro- 
cess starts. It is hence desirable to compute in advance the shortest replacement 
paths, for the case of a failure of any one of the edges in the nominal path P. 
In the event of a failure the routing mechanism simply chooses the correspond- 
ing pre-computed path. The applicability of RP is, however, limited to those 
situations, in which it is possible to know the identity of the failed link before 
the routing process starts. This assumption is not realistic in many important 
applications, in which faults in the network occur 'online', or the information 
about them is stored in a distributed fashion. The latter situation is common- 
place, for example, in transportation networks (e.g. accidents in road networks). 
Furthermore, it is a common feature of very large networks, such as the Internet. 
This paper studies the online replacement path problem (ORP), which captures 
this online failure setting. The most notable difference between RP and ORP is 
that in ORP we assume that the routing mechanism is informed about the failed 
link at the moment it tries to use it. Another important difference is related to 
the nominal path P. In RP a certain nominal shortest path is provided in the 
input. This is no longer the case in ORP, since the detour taken in the event of a 
failure does not always start from the source vertex s. This means that in ORP 
we simultaneously optimize both the nominal path and the optimal detours, tak- 
ing into account a certain global objective function. An informal formulation of 
ORP is as follows. We would like to route a certain package through a network 
from a given source to a given destination as quickly as possible. We are aware 
of the existence of a failed link in the network, but we do not know its location. 
It is possible to observe that a certain link has failed by probing it. In order to 
probe a link, the package should be at one of the endpoints of this link. If a 
probed link is intact, the package crosses the link to the other endpoint and the 
cost of traversing the link in incurred. Otherwise the package stays in the same 
endpoint and the routing mechanism is informed about the failed link. In other 
words, it is only possible to observe that a link has failed by trying to cross it. 
The goal is to find a set of paths (a nominal path and detours for every edge 
on the nominal path) that will minimize the latest possible arrival time to the 
destination. A solution of ORP should hence specify both the nominal path and 
the optimal detours at every vertex along the path, which avoid the next edge 
on the nominal path. 

The latter informal definition suggests that in ORP we take a conservative 
fault-tolerant approach. In fact, it is assumed that a failed link does exist in 
the network, but its identity is unknown. This suggests another application of 
ORP. In many applications it is only necessary to route a certain object within 
a certain time, called a deadline. As long as the object reaches its destination 
before the deadline, no penalty is incurred. On the other hand, if the deadline 
is not met, a large penalty is due. An example of such an application is organ 
transportation for transplants (see e.g. Moreno, Vails and Ribes [2]), in which it 
is critical to deliver a certain organ before the scheduled time for the surgery. In 
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this application it does not matter how early the organ arrives at the destination, 
as long as it arrives in time. In such applications it is often too risky to take 
an unreliable shortest path, which admits only long detours in some scenarios, 
whereas a slightly longer path with reasonably short detours meets the deadline 
in every scenario. Hence, with ORP it is often possible to immunize the path 
against faults in the network. 

The main result of this paper is that the solution to ORP in undirected net- 
works can be computed in 0(mlog?i) time and linear space for all sources and 
a single destination, where n and m are the number of vertices and edges in the 
network, respectively. Furthermore, this solution can be stored in 0(n) space. 
The basic algorithm is a label-setting algorithm, similar to Dijkstra's algorithm 
for ordinary shortest paths. We describe this algorithm in Section [3l The main 
technical difficulty lies in the need to pre-compute certain shortest path distances 
in modified graphs. This difficulty is overcome in Section|4j which provides a fast 
implementation of this step. We generalize the model to incorporate the possi- 
bility of an arbitrary bounded number k of failed edges in Section [5] This section 
also gives an alternative view on the problem using the notion of routing strate- 
gies, and provides a polynomial algorithm for the problem. We mention other 
results linking ORP with the shortest path problem in Section [6l In particular, 
we show that it is possible to solve a bi-objective variant of ORP. This result 
implies that it is possible to efficiently obtain Pareto-optimal paths with respect 
to ordinary distance and the cost corresponding to the ORP problem, making it 
an attractive method for various applications. We also analyze the performance 
of a greedy heuristic which attempts to route along a shortest path in the re- 
maining network. We show that this heuristic, which is implemented in many 
applications, is a poor approximation for ORP. We summarize in Section [7J The 
following section reviews related work. 

2 Related Work 

The replacement path problem was first introduced by Nisan and Ronen [T5] . 
The motivation for their definition stemmed from the following question in auc- 
tion theory: what is the true price of a link in a network, when we try to con- 
nect two distinct vertices x and y, and every edge is owned by a self-interested 
agent. It turns out that compensating agents with respect to the declarations of 
other agents in the auction leads to truthful declarations, namely declarations 
which reflect the true costs of the agent. Such pricing schemes are called Vickrey 
schemes. In the setup of networks, replacement paths lengths correspond to Vick- 
rey prices for the individual failed edges. Another important application of RP 
is the k shortest simple paths problem (kSSP), which reduces to k replacement 
path computations. 

The complexity of the RP problem for undirected graphs is well understood. 
Malik, Mittal and Gupta [T3] give a simple 0(m+n log n) algorithm. A mistake in 
this paper was later corrected by Bar-Noy, Khuller and Schieber [Jj . This running 
time is asymptotically the same as a single source shortest path computation. 
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Nardclli, Proietti and Widmayer [T7] later provided an algorithm with the same 
complexity for the variant of RP, in which vertices are removed instead of edges. 
The same authors give efficient algorithms for finding detour-critical edges for a 
given shortest path in |15I16] . 

In directed graphs the situation is significantly different. A trivial upper 
bound for RP corresponds to 0(n) single shortest path computations. This gives 
0(n(m + nlogn)) for general directed graphs with nonnegative weights. This 
was slightly improved to 0(mn + n 2 log log n) by Gotthilf and Lewenstein [TU] . 
The challenge of improving the O(mn) bound for RP on directed graphs was 
mainly tackled by restricting the class of graphs or by allowing approximate 
solutions. Along the lines of the former approach, algorithms were developed for 
unweighted graphs (Roditty and Zwick [2Tj ) and planar graphs (Emek, Peleg 
and Roditty [7], Klein, Mozes, and Weimann [12] and Wulff-Nilsen [23]). The 
latter approach was successfully applied to obtain |-approximate solutions by 
Roditty [20] and (1 + e)-approximate solutions by Bernstein [5]. Weimann and 
Yuster [22] applied fast matrix multiplication techniques to obtain a randomized 
algorithm with sub-cubic running time for certain ranges of the edge weights. 

Another problem which bears resemblance to ORP is the stochastic shortest 
path with recourse problem (SSPR), studied by Andreatta and Romeo [3]. This 
problem can be seen as the stochastic analogue of ORP. 

Finally, we briefly review some related work on robust counterparts of the 
shortest path problem. The shortest path problem with cost uncertainty was 
studied by Yu and Yang [53], who consider several models for the scenario set. 
These results were later extended by Aissi, Bazgan and Vanderpooten [2] . These 
works also considered a two-stage min-max regret criterion. Dhamdhere, Goyal, 
Ravi and Singh [6] developed the demand-robust model and gave an approxima- 
tion algorithm for the shortest path problem. A two-stage feasibility counterpart 
of the shortest path problem was addressed in Adjiashvili and Zenklusen pQ. 
Puhl [19] provided hardness results for numerous two-stage counterparts the 
shortest path problem, and gave some approximation algorithms. 

3 An Algorithm for ORP 

In this section we develop an algorithm for ORP. Some technical proofs are left 
in Appendix [A] Let us establish some notations first. We are given an undirected 
edge- weighted graph G — (V,E,£), a source s € V and destination t E V. We 
are assuming throughout this paper that the edge weights I are nonnegative. 
For two vertices u,v £ V let V u , v denote the set of simple u-v paths in G. 
Let N(u) denote the set of neighbors of u in G. For a set of edges A C E let 
= J2e£A^( e )- For an edge e G E and a set of edges F C E, let G - e 
and G — F denote the graph obtained by removing the edge e and the edges in 
F, respectively. For a graph H let dn(-, •) denote the shortest path distance in 
H . Paths are always represented as sets of edges, while walks are represented 
as sequences of vertices. For a path P with incident vertices u and v let P[u,v] 
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denote the subpath of P from u to v. For an edge e G P and u G V let 

s u e = d G - e (u,t) 

denote the shortest u-t path distance in G — e. 

Our algorithm uses a label-setting approach, analogous to Dijkstra's algo- 
rithm for shortest paths. In other words, in every iteration the algorithm updates 
certain tentative labels for the vertices of the graph, and fixes a final label to a 
single vertex u. This final label represents the connection cost of u by an optimal 
path to t. 

Definition 1. Given a vertex v G V , the robust length of the v-t path P is 
Val(P) = max{^(P), max {l(P[v,u}) + s^"'}}. 

uu'&P 

The potential y(v) is defined as the minimum of Val(P) over all P G Vvj, and 
any path P* attaining Val(P*) = y(v) is called an optimal nominal path. Finally, 
ORP is to compute y(s) and obtain a corresponding optimal nominal path. 

The robust length of a v-t path P is simply the maximal possible cost incurred 
by following P until a certain vertex, and then taking the best possible detour 
from that vertex to t which avoids the next edge on the path. To avoid confusion, 
we stress that in ORP we assume the existence of at most one failed edge in the 
graph. Consider next a scenario in which an edge uu' G P fails and let u G V be 
the vertex which is closer to v. We can assume without loss of generality that 
the best detour is a shortest u-t path in the graph G — uu'. Critically, the values 
s~ uu are independent of the chosen path. 

Observe that from non- negativity of I we obtain Val(P) > Val(P'), whenever 
P and P' are u-t and v-t paths respectively, and P' is a subpath of P. We denote 
this property by monotonicity. Furthermore, we can prove the following. 

Lemma 1. Let P u G V u ,t and let v G N(u) be a vertex, not incident to P u . 
Then the path P v = P u U {vu} satisfies 

Val(P„) = m&x{e(vu) + Val(P„), s~ vu }. 

Our algorithm for ORP updates the potential on the vertices of the graph, using 
the property established by the following lemma. 

Lemma 2. Let U G V, with t G U, be the set of vertices for which the potential 
is known. And let uv be the edge such that: 

uv — argmin {max{l(zw) + y(w), s~ zw }}. (I) 

zwt£E:weV ,z€V\U 

Then if P u G V u t with Val(P u ) = y(u) and P v = P u U {uv} it holds that 
Val(P„) - y(v). 
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Algorithm 1 

1: Compute s~ uv for each uv 6 E. 

2-.U = 9; W = V; y'(t) = 0; y'{u) = oo \/u e V - 1. 
3: successor(u) — NIL Vu £ V. 
4: while U # V do 
5: Find u = arg min z6W , y'(z). 
6: f/ = t/ + u; W = W-u; y(u) = y'(it). 
7: for all vu £ E with u £ do 
8: if y'(u) >max{fH + i/'(ii),s;™} then 

9: ?/0) = max{f(uu) + y'(u), s~ vv }. 

10: successor(v) — u. 



Lemma [2] provides the required equation for our label-setting algorithm, whose 
formal statement is given as Algorithm [TJ The algorithm iteratively builds up a 
set U, consisting of all vertices, for which the correct potential value of y(u) was 
already computed. The correctness of the algorithm is a direct consequence of 
Lemma [2] 

Consider the running time of Algorithm^ We let n and m denote the number 
of vertices and edges of the input graph, respectively. An efficient implementation 
of step 1 is delayed to the next section, and constitutes the heart of our efficient 
algorithms. For steps 2-10 we use the implementation of Fredman and Tarjan [5] 
for priority queues (heaps) called Fibonacci Heaps. We adopt here the same 
implementation that is used to obtain 0(m + rologn) running time for Dijkstra's 
algorithm. We omit the details as they are identical to those in [8] . We comment 
that our implementation of step 1 uses 0(m log n) time, hence this computation 
dominates the running time. In fact, the running time of Algorithm [T] remains 
the same if steps 2-10 are implemented using simpler data structures, such as 
Binary Heaps. 

Finally, let us remark that Algorithm [T] works both for directed and undi- 
rected graphs. However, the following section provides an implementation of 
step 1 for undirected graphs only. We remark that in directed graphs, Algo- 
rithm[T]can be trivially implemented in time 0(nSP(n, m)), where SP(n, m) is 
the complexity of a single shortest path computation on a graph with n ver- 
tices and m edges. This is a simple consequence of our first observation in the 
following section. 

4 Efficient Computation of s- Values 

It remains to provide an efficient implementation for the computation of the val- 
ues s~ uu . We will assume a random access machine (RAM) as a computational 
model. 

We start by computing the shortest path tree T in O(mlogn) time from 
every vertex to t. Let d*(u) = dc(u,t). Observe that s~ uu = d*(u) holds for 
every edge uu 1 outside of T. We can hence concentrate our efforts on computing 
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s- values for edges in the tree. For a vertex v G V we denote by T v C T the 
subtree rooted at v. We set E' = E\T. 

Lemma 3. Consider a vertex u £ V . Let vw G E' be an edge attaining 

min ltfc(t*,«) +£(vw) +d*(w)}, (2) 

vw£E' ,v£T u w£T lt 

and let uu' be the first edge on the u-t path in T . Then 

s~ uu ' = d G (u, v) + l{vw) + d* (to). 

By defining c vw = d*(v) + d* (w) + £(vw) for each edge vw G E' and substituting 
in the expression for s~ uu obtained in the previous lemma we can write 

C"' =c vw -d*{u), (3) 

which is a convenient expression for computing the s-values. Indeed this expres- 
sion shows that value s~ uu only depends on the lowest value c e , over all e with 
exactly one incident vertex in T u . 

We leave the remaining details of our implementation to the proof of Theo- 
rem [1] Before stating the theorem, let us define an important ingredient, which 
is used hereafter. To efficiently obtain the information, whether a certain edge 
e G E' corresponds to a feasible detour for u (namely if exactly one endpoint of e 
is in T u ), we need an algorithm for computing the least common ancestor (LCA) 
in trees. Given a tree rooted at a vertex t and two vertices u, v in the tree, the 
least common ancestor lca(u, v) of u and v is the vertex at which the t-u and t-v 
paths diverge in the tree. For an edge e we write lca(e) to denote the LCA of 
the endpoints of e. It is straightforward to see that an edge e = vw corresponds 
to a feasible detour for an edge uu' if and only if either the v-\ca,(v,w) path or 
the w-lca(v,w) path in T contain uu' . For this purpose we use the algorithm 
of Gabow and Tarjan [5], which uses 0(n + p) time to compute the LCA of p 
pairs of vertices in a rooted tree. In our case we have p = O(m), since we would 
like to compute the LCA for every pair of vertices connected by an edge e G E', 
hence this computation takes 0(m) time. We assume henceforth that given an 
edge e G E' we have access to lca(e) in constant time. 

Theorem 1. The values s~ uu can be computed in 0{m log n) time and linear 
space for all uu' G E. Furthermore, they can be stored in a data structure of size 
O(n). 

Proof. The algorithm is summarized as Algorithm [5] in Appendix [Bj The algo- 
rithm starts by sorting the set of edges E' according to increasing order of c e . 
Let L be the sorted list. This operation (as well as the computation of the values 
c e ) takes 0(m log n) time. 

The algorithm relies on the following fact. Let e be the first edge in the 
sorted list L, which represents a feasible detour for uu' G T. From the previous 
discussion, the edge e corresponds to the optimal detour for uu' . Furthermore, if 
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we traverse the two paths from the endpoints of e to lca(e) we cross the edge uu' . 
The algorithm indeed performs the latter traversals in a copy of T and marks 
all edges along the way (which were not marked before, using an edge that has 
lower c- value) as belonging to the edge e. To obtain the desired running time we 
need to avoid traversing the same parts of T several times. To achieve this we 
create a second copy T of the tree T. For every vertex u e T we store a pointer 
p[u] from it e T to its copy u in T, and another pointer p[u], pointing in the 
opposite direction. We start iterating over the edges in the sorted list L. In the 
beginning of the i'th iteration, the first edge ej = xy in L is removed from L, and 
the algorithm jumps to the corresponding vertices x = p[x] and y — p[y] in T. 
The algorithm also memorizes z = lca(5, y) - the vertex in T, corresponding to 
z = lca(x, y). Next the x-z and y-z paths are traversed in T, marking all edges 
along the way as belonging to ej, and computing the s-values for them using 

. Finally all these aforementioned edges are contracted in T and the pointers 
are updated so that every vertex in T always points to its corresponding super- 
vertex in T and vice-versa. Finally, the list is post-processed to remove some 
unnecessary edges in the following way. At every step, the first edge e = vw in 
the list is inspected. If p[v] = p[w\, this edges is a self- loop in T, and it can not 
correspond to correct s-value assignments anymore. Consequently, this edges is 
removed from the list and the next edge is inspected. The process ends when 
p[v] ^ p[w], or when L is empty. In the latter case the algorithm terminates. 
This concludes the i'th iteration. 

Note that the running time of the i'th iteration is proportional to the com- 
bined length of the x-z path, the y-z path and the prefix of L that is deleted in 
the post-processing. Since we contract every edge we traverse, in the tree T, we 
never traverse the same edge twice. We conclude that the running time of this 
last stage of the algorithm is 0(n + to). Figure Q] illustrates the algorithm. 

Finally, note that we can store a representation of all s-values and the corre- 
sponding detours in 0(n) space. This is achieved by storing the tree T, and for 
every vertex u and the corresponding edge uu' in T, we store a pointer to the 
edge attaining the minimum in ([2]). The values s~ uu can either be stored in a 
separate list using O(to) additional space, or computed in constant time from 
the aforementioned information, by using (|3|). 

□ 

Note that the running time of Algorithm [5] is dominated by the sorting of the 
costs of the edges. This fact is a significant advantage since sorting algorithms 
are very efficient in practice. In fact, if the sorted list L was provided in the 
input, the running time of the algorithm could be improved to 0(m + nlogra), 
via the implementation mentioned in Section [3J Another case in which the latter 
complexity bound can be attained is that of unweighted graphs, where bucketing 
can be used in order to sort the costs c e in linear time. Theorem [2] summarizes 
our main result. 

Theorem 2. Given an instance of ORP the potential y and the corresponding 
paths can be computed in time 0(m log n). 
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Fig. 1. Top left: The shortest path tree T. The next edge in E' to be processed is 
ej = xy. The cost c xy associated with this edge is the sum of weights of all dashed 
edges, with the edge vt counted twice. lca(xy) = v. Top right: The tree T in the 
beginning of the i'th iteration. The labels of the super-vertices correspond to the sets 
vertices in T associated with them. The dashed edges are traversed in the i'th iteration 
and contracted to form the tree on the bottom of the figure. 

5 fc-ORP 

Let us formally define fc-ORP, the online replacement path problem with k failed 
edges. We refer to the parameter k as the failure parameter. In /c-ORP a scenario 
corresponds to a removal of any k of the edges in the graph. In this setup it is no 
longer convenient to describe the problem in terms of paths and detours. Instead 
we introduce the notion of a routing strategy. 

A routing strategy R : 2 E x V — > V is a function which, given a subset 
F' C E of known failed links and a vertex v £ V, returns a vertex u G V. We 
call E' = E\F' the set of active edges. We are assuming the existence of a certain 
governing mechanism, which takes as input a routing strategy and executes it on 
a given instance. Provided with a routing strategy R, this mechanism iteratively 
moves from the current vertex u to the vertex v = R(F', u), where F' is the set 
of failed links probed so far. The process starts at a given origin s with F' = 
and ends when t is reached. 

Since R is deterministic, this process defines a unique, possibly infinite, walk 
6r(u, E, F) in G for each origin u S V and every scenario F C E with \F\ < k. 
We remark that if G contains at least k + 1 edge-disjoint s-t paths, there exists 
a routing strategy, which does not cycle. 
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Definition 2. Given E' C E, a vertex u £ V and a routing strategy R, the 
k-value of R with respect to E' and u is defined as 

Vfdk(u,E',R)= max £(9 R (u, E' , F)). 

FQE',\F\<k 

The corresponding k-potential of u with respect to E' is defined as 
y k (u, E') = min Val fe (u, E', R). 

Finally, fc-ORP is to find an optimal routing strategy R*, which minimizes 
Vah(s, E, R), namely to solve 

R* = argmin Vah(s, E, R). 

R 

The following relation, which is a generalization of ([1]), gives rise to a simple 
recursive algorithm for fc-ORP. For each v 6 V it holds: 

y k (v,E)= min {ma,x{£(vu) + y k (u,E) 1 y k ~ 1 (v,E\vu)}}. (4) 

v:vu£E 

The correctness of this relation can be proved by induction on k, using the argu- 
ments in Section[3] The exact statement of the algorithm is given in Appendix [Cl 
alongside a complexity analysis. We summarize with the following theorem. 

Theorem 3. k-ORP can be solved on undirected graphs in time 0(m k \ogn). 
6 ORP vs. Shortest Paths 

As we have seen, the length of a path and its robust length are different, and 
often conflicting objectives functions. In this section we mention two results 
relating these two objectives. 

Consider first the problem of finding an optimal solution to an instance of 
ORP with the shortest possible nominal path length. This bi-objective problem 
asks to find a Pareto-optimal path with respect to the latter two objective func- 
tions. Our first result asserts that this problem can be solved in polynomial time 
with a simple adaptation of Dijkstra's algorithm. In fact, we are able to solve 
the following problem for every bound B > OPT. 

P* = argmin l(P)- 

PeV s ,f- Val(P)<B 

The algorithm and possible applications of this problem are given in Appendix iDl 
The second result analyzes the performance of certain greedy heuristics for 
fc-ORP. We show that a routing strategy that always tries to route along the 
shortest path in the remaining graph is a (2 k+1 — l)-approximation algorithm 
for fc-ORP. The details are available in Appendix [El 
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7 Conclusions 

This paper introduces a natural variant of the replacement path problem, the 
online replacement path problem. ORP captures many real-life situations, which 
occur in faulty or large distributed networks. The most important characteristic 
of ORP is that the information about the failed edge in the network is only 
available locally, namely on the edge itself. ORP possesses many of the nice 
characteristics of ordinary shortest paths. In particular, we show that ORP can 
be solved by a simple label-setting algorithm. The computational bottleneck of 
this algorithm is the need to pre-computc certain shortest paths in an adapted 
network. We give an efficient implementations of this step, which uses 0(m log n) 
time and linear space. We also generalize the algorithm to deal with an arbitrary 
constant number k of failed edges. Along this vein we introduce the notion of a 
routing strategy. Finally, we observe that a Pareto-optimal path with respect to 
ordinary distance and robust length can be found in polynomial time. 

We conclude by mentioning a number of promising directions for future re- 
search. The complexity of ORP with a variable number k of failed edges remains 
open. In fact, it is not clear if the decision problem y k (s,E) < M is in NP. 
The complexity of ORP in undirected graphs may potentially be improved to 
0(77j + nlogn). Finally, the complexity of ORP in directed graphs remains open. 
Our results give an algorithm with the running time of 0(nSP(n, m)), where 
SP(n, m) is the complexity of a single shortest path computation on a graph 
with n vertices and m edges. This problem threatens to be as challenging as RP 
in directed graphs. 
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A Proofs 

A.l Proof of Lemma [1] 

Applying the definition we compute 

t(vu) + Val(P„) = i(vu) + max{j?(P„), max„ 6P „ {£(P u [u, z]) + s- zw }} 

= max{!(ra) + £[P U ), max zto£P J!(w) + £{P u [u, z}) + s~ zw }} 
= max{e(P v ),max zwePv \ {vu} {f(P v [v,z]) + s~ zw }}. 

Substituting in the desired expression we obtain 
max{!(i)ii) + Val(P u ), s~ vu } = 

ma X {e{P v ),max zwePA{vu} {e(P v [u,z]) + sJ zw },s- vu } = Val(P„), 
which proves the lemma. 

A. 2 Proof of Lemma [2] 

Assume towards contradiction that there exists a path P* G V Vj t such that 
Val(P*) < Val(P„). Let w EU and z e V \ U be two vertices such that zw G P* . 
Consider the partition of P* given by P„* = P*[v, z] U {zai}UP„*[w), t] (Note that 
if io = t then P*[iy, t] = and if z = v then P*[v, z] = 0). 
By the choice of we have 

max{f( TO ) + y(u), s'™} < max{f(z ffl ) + y(w), S - zw }, 

which, by Val(P u ) = y(u) and Val(P*[w, i\) > y(w) implies 

max{^(wu) + Val(P u ), S -™} < max{f(zw) + Val(P; [w,t]), s~ zw }. 

The latter inequality and Lemma [1] give Val(P„) < Val(P* [z , i\) . On the 
other hand, by monotonicity we have Val(P*) > Val(P* [z, t}). We conclude that 
Val(P„) < Val(P*[z,t}) < Val(P*); a contradiction. 

A. 3 Proof of Lemma [3] 

Clearly this choice represents a u-t path in G—uu', so s~ uu < dc(u 7 v) +£(vw) + 
d*(w) holds. Assume towards contradiction that a better u-t path P existed in 
G — uu' , namely £(P) < da(u,v) + £(vw) + d*(w). To reach t from u without 
using uu 1 the path P needs to contain an edge v'w', such that v' G T u and 
w' <£T U . We obtain 

1{P) =£(P[u, v'])+ £{v'w')+l{P[w',t}) > d G [u,v'}+e(v'w')+d*(w') 

> d G (u,v) +£(vw) + d*(w), 

which contradicts the choice of P. 
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B A Summary of the Algorithm in Section HI 



Algorithm 2 

1: Compute the shortest path tree T. 
2: Create a copy T of T. 

3: Store a pointer p[u] in every vertex u in T to its corresponding vertex u in T. 
4: Store a pointer p[u] in every vertex u in T to its corresponding vertex u in T. 
5: Store similar pointers p[e] from edges in T to corresponding edges in T. 
6: Compute c e for every e <E E' . 

7: Sort iJ' according to increasing order of c e . Let L be the sorted list. 

8: while T contains more than one vertex do 

9: Remove the first edge e = xy from L. 
10: x = p[x], y=p[y]. 
11: z = lca(e), z = p[lca(e)]. 

12: Let Y <ZT denote the union of the x-z and y-z paths in T. 

13: for a G Y do 

14: uu — p[a], 

15: s"""' =Ce-d*(u). 

16: Contract F in T and update pointers. 

17: repeat 

18: Let e = vw be the first edge in L. 
19: if p[«] = p[w] then 

20: Remove e from L. 

21: until p[v] ^ p[w] or L = 0. 



C An Algorithm for fc-ORP 

Algorithm [3] solves (J4j) for every v E V ~ t and computes the corresponding 
optimal routing strategy. Note that Algorithm [3] is identical to Algorithm Q] 
when k = 1 . To formally prove the correctness of Algorithm [3] one needs to 
state equivalent monotonicity properties, as well as analogues of Lemma [T] and 
Lemma [2j They are however omitted, as they are identical to those of Section [3l 
We conclude by bounding the complexity of a naive implementation of 
this algorithm. Let T(m,n,k) denote the running time of the algorithm on a 
graph with n vertices and m edges and failure parameter k. The algorithm 
in Section U gives T(n,m, 1) = 0(m log n). For k > 1 we have T(m,n,k) = 
0(m + n logn) + mT(m — 1, n, k — 1) = 0(m k ~ 1 T(n, m, 1)). This finally gives 
T(n,m, k) = 0(m k logn). 

D Pareto-Optimal ORPs 

In this section we are concerned with obtaining a path P* with robust length at 
most B, and a shortest nominal path length among all such paths. This problem 
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Algorithm 3 



1: Compute y k 1 (u, E\ uv) for each uv 6 E. 

2: (7 = 0. 

3: W = V. 

4: y k (t,E) = 0. 

5: y k (u,E) = oo for each u £ VF. 

6: while U / V do 

7: Find it = argmin zSVK {y fc (>, £)}. 

8: U = U + u. 

9: W = W-u. 

10: for all vu £ E such that n € W 7 do 

11: if y k (v,E) > ma,x{e(uv) + y k (u,E),y k - 1 (v,E\{vu})} then 

12: j^O, S) = max{£(™) + y h (u, E), y^iv, E \ {vu})}. 



is equivalent to that of finding a Pareto-optimal path with respect to the objec- 
tive functions corresponding to the ordinary distance function, and the robust 
length. We assume B > OPT, the optimal solution value of the corresponding 
ORP instance. Formally, we aim at finding an s-t path P* satisfying 

P* = argmin l(P)- 

PeV,,f. Val(P)<B 

The algorithm for this problem is a slightly modified version of Dijkstra's 
algorithm on the directed edge- weighted graph G 1 = (V,K,lu), where K has 
two directed edges uv and vu for each undirected edge uv £ E, with ui(uv) = 
ui(vu) = £(uv). Algorithm |4] is a formal statement of the algorithm. Note that 
the only difference with Dijkstra's algorithm is the condition in step 8. Unlike 
ordinary shortest paths, when a vertex u is selected, and the distance labels of 
its neighbors are updated, it is not sufficient to check for each uv G K the usual 
condition 

d(u) + oj(uv) < d(v). 

We need to additionally verify whether the edge uv can belong to a path with 
robust length B or not. In other words, we need to check whether the length of 
the path from s to u plus the length of the detour from u to t avoiding uv is at 
most B, namely 

d(u) + s~ uv < B. 

Algorithmic can clearly be implemented to run in 0{m + n log n) time using the 
results of Section |U 

Let us briefly discuss the potential applications of the latter problem. While 
shortest paths often have undesirable behavior in unreliable networks, an optimal 
solution to ORP might have a prohibitively large cost. In some applications faults 
occur rarely, hence it is preferred to have the cost of the nominal path as low 
as possible. At the same time, it is necessary that the cost does not exceed a 
certain threshold B, in every scenario. Consequently, it makes sense to regard 
the threshold B as a hard constraint on the robust length, and optimize the 
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Algorithm 4 

1: S = <t); S = V 

2: d(s) = 0; = oo \/v G V - s 

3: while £ £ 5 do 

4: Find u = argmin ze 5 

5: S = S + u 

6: S = S — u 

7: for w G iV(u) \Sdo 

8: if d(u) + uj{uv) < d(v) and d(u) + s~ wv < B then 

9: d(v) = d(u) + u>(uv) 



length of the nominal path. Algorithm @] gives the decision maker the desired 
freedom to choose the level of conservatism that she desires. 

E The Shortest Path Heuristics 

It is common to use heuristics which rely on shortest paths in routing algorithms. 
In this section we show that a naive shortest path routing strategy performs 
very poorly for fc-ORP. In fact, the approximation guarantee it provides grows 
exponentially with the adversarial budget k. To this end we define more formally 
the shortest path heuristics, which we denote by Rsp- The routing strategy 
Rgp works as follows. At each vertex u G V and given a set of known failed 
edges F', Rsp tries to route the package along a shortest path in the remaining 
graph G — F' . In the following lemma we show that Rsp is a factor 2 fc+1 — 1 
approximation for the optimal routing strategy, in the presence of at most k 
failed edges. 

Lemma 4. Let I = (G,s,t) be an instance of k-ORP. Then 

Val fc (s, E,R SP ) < (2 fc+1 - l)y k {s,E). 

Proof. Let OPT = y k (s,E), R* = arg min fl Vah (s, E, R) and Q° = 9r.(s,EJ) 
be the corresponding nominal path. Consider a set F of failed edges with |F| < k. 
Define ui = 8r sp (s, E, F) to be the walk followed by Rsp in G — F, and let 
lu = (s = u±,U2, ■■■ ,u m = t) be the corresponding sequence of vertices. We 
divide our analysis according to the number of failed edges encountered by Rsp- 
We prove by induction on i, that if the routing strategy encountered a total of 
i < k failed edges then 

t(u) < (2 i+1 - l)OPT. (5) 

The base case i = corresponds to scenarios in which no failed edge is encoun- 
tered in the routing. In this case oj is simply a shortest s-t path in G, hence 
l{uS) < OPT, as required. Assume next that ([5]) holds for every j < i and con- 
sider the case that the routing encounters exactly i failed edges. We can assume 
without loss of generality that all edges of F were probed and |F| = i. 
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Let r < m be such that u r is the vertex incident to the i'th failed edge 
in the routing. In other words, before reaching u r , the routing probed exactly 
i — 1 failed edges. Let e = {u r , w} be the i'th failed edge probed by the routing 
strategy. Consider an execution of the routing strategy on the same instance 
with the different failure scenario F 1 = F — e. The resulting walk uj' will have 
the first r vertices in common with u>, namely the sub- walk a = (s = u\, ■ ■ ■ ,u r ) 
will appear in both walks. By the inductive hypothesis we have that £(a) < 
^(w') < (2* — \)OPT. It remains to bound the length of the tail of ui from u r 
until u m = t to complete the proof. To this end recall that Rsp routes the 
package along the shortest remaining path in the graph. Since the last failure 
encountered by Rsp is e, the remaining path is simply the shortest u r -t path 
in G — F. To bound the length of this path we construct a u r -t walk 9 as 
follows. First 9 traces the entire route taken by Rsp back to s and then uses 
the walk that R* would use to reach t from s in the scenario F. Clearly we have 
l{6) < (2 l - l)OPT + OPT. Furthermore, this walk is intact in G-F. This gives 
the required bound < (2 1 - l)OPT +(T — l)OPT + OPT = (2 l+1 -l)OPT 
and finishes the proof. □ 

The bound obtained in Lemma|4]seems crude at first glance. In particular, in 
the inductive step we follow the entire walk performed so far backwards to reach 
s and start over. In the following example we show that the bound of Lemma 0] 
is tight. 

Example 1. Let M G Z + be a large integer. Consider the following instance 
I = (G, s, t) of fc-ORP. The graph contains k + 1 parallel edges connecting s and 
t with length M + 1. In addition the graph contains a path (s,ui, ■ ■ ■ ,Uk) of 
length k + 1. The edge sui has length M and every edge UiUi + i has length 2 l M. 
Finally, the vertices u\, ■ ■ ■ ,Uk are connected to t with edges of length zero. The 
construction is illustrated in Figure [2j 

Consider the failure scenario, which fails all edges u{t for i £ [k]. The routing 
strategy Rsp will follow the path (s, Ui, • ■ • then follow it back to s and 

then take one of the edges with length M + 1 to t. The total length of this walk 

is 2(M + 2M H h 2 k M) + M + 1 = (2 fc+1 - 1)M + 1. At the same time the 

optimal routing strategy routs the package along the edges with length M + 1 
with a worst-case cost of M + 1. The ratio between the two numbers tends to 
2 fe+1 — 1 as M tends to infinity. 
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